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Abstract 

It is shown that the capture and fission cross-sections of 238 U and 239 Pu increase 
with temperature within 1000-3000 K range, in contrast to those of 235 U, that under 
certain conditions may lead to the so-called blow-up modes, stimulating the anomalous 
neutron flux and nuclear fuel temperature growth. Some features of the blow-up regimes 
in neutron-multiplying media are discussed. 

1 Introduction 

It is known that with loss of coolant at three nuclear reactors during Fukushima nuclear 
accident its nuclear fuel melted. It means that temperature in the active zone in some 
moments reached melting point of uranium-oxide fuel 1 , i.e. ^3000°C. 

Surprisingly in scientific literature today there are absolutely no neither experimental 
nor even theoretically calculated data on behavior of capture and fission cross-sections of 
238 U and 239 Pu depending on temperature at least in 1000-3000°Crange. At the same 
time there are serious reasons to suppose that the cross-section values of the specified 
elements increase with temperature increase. At least we may point, for example, to the 
qualitative estimates by Ukraintsev [1], Obninsk Institute of Atomic Energetics (Russia), 
that confirm the possibility for cross-sections growth phenomenon for 239 Pu in 300-1 500° C 
range. 

Obviously, such anomalous temperature dependence of capture and fission cross-sections 
of 238 £7 and 239 Pu may drastically change neutron and thermal kinetics of nuclear reactors 
including perspective fast new generation reactors (uranium-plutonium of Feoktistov [2] 
and thorium- uranium of Teller [3] type), that we classify as fast TWR reactors. Hence 
it is very important not only to know the anomalous temperature behavior of 238 U and 
239 Pu capture and fission cross-sections, but also it becomes critically important to know 
their influence on heat transfer kinetics, because it may turn into a reason of positive 
feedback 2 (PF) with neutron kinetics leading to undesirable solution stability loss (the 

* Corresponding author: Vitaliy D. Rusov, E-mail: siiis@te.net.ua 

1 We should note here, that the third block partially used MOX-fuel, enriched with plutonium 

2 Positive Feedback is a type of feedback when change in output signal leads to the change in input signal, 

that in its turn leads to further deviation of output signal from its original value. In other words PF leads to 

instability and appearance of qualitatively new (often self-oscillation) systems. 
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nuclear burning wave) as well as to a trivial reactor runaway with a subsequent nontrivial 
catastrophe. 

A special case of PF is a non-linear PF, leading to system evolution in so-called blow- 
up mode [4-9] or, in other words, in such a dynamic mode when one or a number of 
modeled values (e.g. temperature and neutron flux) grows to infinity at a finite time. In 
reality instead of infinite values a phase transformation is observed in this case, that can 
in its turn become a first stage or a precursor of future technogenic disaster. 

Investigation of temperature dependence of capture and fission cross-sections of 238 U 
and 239 Pu in 300-3000° C range and, correspondingly, the kinetics of heat transfer and its 
influence on neutron kinetics in TWR is the main goal of the present paper. 

2 Temperature blow-up regimes in neutron-multiplying 
media 

Heat transfer equation for uranium-plutonium fission environment is: 

dT (r, t) 



p{r,T,t).c{r,T,t) 



dt 



= N (f, T, t) • AT (f, t) + W (f, T, t) • VT (r, t) + q J T (f, T, t) , (1) 
where the effective substance density is 

P (r,T,t) = Y^Ni(r,T,t)- P i, (2) 

i 

Pi are tabulated values, N{ (r, T, t) are the components concentrations in the medium, 
while the effective specific heat capacity ( accounting for the medium components heat 
capacity values q) and fissile material heat conductivity coefficient (accounting for the 
medium components heat conductivity coefficients respectively are: 

c(f,T,t) = Y,c i (T)N i (f,T,t), (3) 

i 

X(r,T,t) = J2HT)Ni(r,T,t). (4) 

i 

Here (r , T, t) is the heat source density generated by nuclear fissions N{ of fissile 
metal composition components that vary in time. 

Theoretical temperature dependency of heat capacity c(T) for metals is known: at low 
temperatures c(t) ~ T 3 , and at high temperatures c(T) — >► const, notably the constant 
value (const w 6 Cal/(mol • deg)) is determined by Dulong-Petit law. At the same time it 
is known that thermal expansion coefficient is small for metals, therefore the specific heat 
capacity at constant volume c v almost equals to the specific heat capacity at constant 
pressure c p . On the other hand, the theoretical dependency of heat conductivity Ni(T) at 
high temperature of "fissile" metals is not known, while it is experimentally determined 
that heat conductivity coefficient N(T) of fissile medium is a non-linear function of tem- 
perature (e.g. see [10], where heat conductivity coefficient is given for a-uranium 238 and 
for metallic plutonium 239, and also [11]). 

Further for solving thermal conductivity equation we used the following initial and 
boundary conditions: 

T(r, t = 0) = 300 K and j n = « [T(r G », t) - T ] , (5) 
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where j n is the normal (to the fissile medium boundary) heat flux density component, H(T) 
is the thermal conductivity coefficient, $ft is the fissile medium boundary, Tq - temperature 
of the medium adjacent to the active zone. 

Obviously, if the cross-sections of some fissile nuclides increase with temperature, then 
due to exothermic nuclei fission reaction the direct consequence of significantly non-linear 
kinetics of mother and daughter nuclides in the nuclear reactor would be an autocatalytic 
increase of generated heat, similar to autocatalytic processes of exothermic chemical re- 
actions. In this case thermal flux density c^ T (r, $,T, t) that characterizes generated heat 
amount will be: 

4 (f, T, t) = $ (r, T, t) J2 Q{ 4 (r, T, t) N* (f, T, t) , [W/cm 3 ] , (6) 

i 

where 

4>(r,T,t) = J <S>(r,E,T,t)dE 
o 

is the full neutron flux density; $ (r, E,T, t) - neutron flux density with energy E\ Q{ - 
mean generated heat emitted due to fission of one nucleus of the z-th nuclide; 

<f)(r,T,f) = J a}(E,T)p(r,E,T,t)dE 
o 

is the fission cross-section of the z-th nuclide averaged over the neutron spectrum; 

p (f, E, T,t) = $ (r, E, T, t) /$ (f, T, t) 

is the neutron energy distribution probability density function; c l f(E, T) is the microscopic 
fission cross-section of the i-th nuclide that, as known, depends on neutron energy and 
fissile medium temperature (Doppler effect [12]); N{ (r, T, t) - density of the i-th nuclide 
nuclei. 

As follows from (6), in order to build thermal source density function it is necessary to 
solve a problem of theoretical dependence of the cross-sections (r, T, t) averaged over 
the neutron spectrum on the temperature of reactor fuel (fissile medium). As is known, 
the influence of the nuclei thermal motion in the medium comes to a broadening and 
height reduction of the resonances. By optical analogy this phenomenon is referred to 
as Doppler effect [12]. Since in low energy area the resonance levels are observed only 
for heavy nuclei, then Doppler effect is noticeable only during the interaction of neutrons 
with such nuclei. And the higher environment temperature the stronger is the effect. 

Therefore a program was developed using Microsoft Fortran Power Station 4.0 (MFPS 
4.0) that allows at the first stage to calculate the cross-sections of the resonance neutron 
reactions depending on neutron energy taking into account the Doppler effect. The cross- 
sections dependence on neutron energy for reactor nuclides from ENDF/B-VII database [13], 
corresponding to 300K environment temperature, were taken as the input data for the 
calculations. For example, the results for radioactive neutron capture cross-sections de- 
pendence on neutron energy for 235 U are given in fig. 1 for different temperatures of the 
fissile medium in 300K-3000K temperature range. Using this program, the dependence 
of scattering, fission and radioactive neutron capture cross-sections for essential reactor 
fuel nuclides g 35 ^, i! 8 ^? 92 9 ^ an d q^9Pu and for different temperatures in range 300K 
to 3000K were obtained. 

At the second stage a program was developed to calculate the dependence of the 
cross-sections (r, T, t) averaged over the neutron spectrum for main reactor nuclides and 
for main reactor nuclear neutron reactions for the specified temperatures. The averaging 
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Figure 1: Calculated dependency of radioactive neutron capture reaction on its energy for H U 
at different temperatures within 300K to 3000K. 



of neutron cross-sections for the Maxwell distribution was performed using the following 
expression: 

(a (Eum, T)) = — — , 

jf E 1 / 2 e~ E / kT dE 
o 

where En m is the the upper limit of neutrons thermalization, while for the procedure of 
neutron cross-sections averaging over Fermi spectrum the following expression was used: 

oo 

/ a(E,T)E~ 1 dE 

(a(E lim ,T)) = , 

/ E~ x dE 

During further calculations in our programs we used the results obtained at the first 
stage i.e. the dependence of reaction cross-sections on neutron energy and medium tem- 
perature (Doppler effect). The neutron spectrum was specified in a combined way - below 
the limit of thermalization En m the neutron spectrum was described by Maxwell spectrum 
<E>M (En); above Eu m but below Ep (upper limit for Fermi neutron energy spectrum) the 
neutron spectrum was described by Fermi spectrum <&f(E) for a moderating medium with 
absorption; above Ef, but below maximal neutron energy £™ x the spectrum was de- 
scribed by 239 Pu fission spectrum [14, 15]. Here the neutron gas temperature for Maxwell 
distribution was given by (7), described in [12]. According to this approach [12], the draw- 
backs of standard slowing-down theory for thermalization area may be formally reduced 
if a variable = £(1 — 2/z) is introduced instead of the average logarithmic energy loss 
£, which is almost independent of neutron energy (as is known, for environment consisting 
of nuclei with A > 10 the statement £ « 2/ A is true). Here z = E n /kT, E n is the neutron 
energy, T - environment temperature. Then within such a framework for neutron gas 
temperature in Maxwell spectrum of thermal neutrons the following expression may be 
used 1 : 

A very interesting expression revealing hidden connection between temperature of neutron gas and envi- 
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(7) 

where Tq is the fuel environment temperature, S a (/cTo) is an absorption cross-section 
for energy fcTo, 77 = 1.8 — dimensionless constant, (£) is averaged over the whole energy 
interval of Maxwell spectrum £(z) at kT = 1 eV. 

Fermi neutron spectrum for a moderating medium with absorption (we considered 
carbon as a moderator and 238 U, 239 U and 239 Pu as absorbers) was set in the form 
[12, 16]: 

s 

* 'Fermi (E, Ep) = ^^ t E ^ 

where S is the total volume neutron generation rate, (£) = ^ fe^s) /^Sj & i s the 

average logarithmic decrement of energy loss, Sg - macroscopic scattering cross-section 
of z-th nuclide, St = Yl + ~~ total macroscopic cross-section of the fissile material, 

i 

S^ = Yl ~ total macroscopic scattering cross-section of the fissile material, S a - 

i 

macroscopic absorption cross-section, Ep - upper neutron energy for Fermi spectrum. 

The upper limit of neutron thermalization Eu m in our calculation was considered as 
a free parameter, setting the neutron fluxes of Maxwell and Fermi spectra at a common 
energy limit Eu m equal: 

$ Maxwell (Elim) ^ Fermi (Eu m ) . (9) 

The high energy neutron spectrum part (E > Ep) was defined by fission spectrum 
[16-18] in our calculations. Therefore for the total volume neutron generation rate 5, in 
the Fermi spectrum expression (8), the following expression may be written : 



T 



1 + 7? 



l 



En 



S Q (g) dE' 
(C)S t {E') E' 



(8) 



5(r,T,t) 



PftE,T,t) 



E F 



Ui(E) • $ (f, E, T, t) • g){E, T) • Ni (f, T, t) 



dE, (10) 



where E™ ax is the maximum energy of neutron fission spectrum (usually taken as E™ ax « 
10 MeV), Ep is the neutron energy, below which the moderating neutrons spectrum is 
described as Fermi spectrum (usually taken as Ep ~ 0.2 MeV); P (r, E, T, t) - probability 
of neutron not leaving the boundaries of fissile medium which depends on fissile material 
geometry and conditions at its border as well (e.g. presence of a reflector). 

The obtained calculation results show that cross-sections averaged over the spectrum 
may increase (fig. 2 for 239 Pu and fig. 4 for 238 U) as well as decrease (fig. 3 for 235 U) 
with fissile medium temperature increase . As follows from the obtained results, the 
arbitrariness in selection of the limit energy for joining Maxwell and Fermi spectra does 
not significantly alter the character of these dependencies evolution. 

This can be justified by the fact that 239 Pu resonance area starts from significantly 
lower energies than that of 235 JJ and with fuel temperature increase, the neutron gas 
temperature increases producing Maxwell's neutron distribution maximum shift to higher 
neutron energies. In other words, the neutron gas spectrum hardening, when more neu- 
trons fit into resonance area of 239 Pu, is the cause of the averaged cross-sections growth. 

For 235 U this process in not as significant because its resonance area is located at higher 
energies. As a result, the 235 U neutron gas spectrum hardening related to fuel temperature 
increase (in the considered interval) does not result in significant increase of a number 
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Figure 2: Temperature dependencies for fission cross-section (a) and radioactive capture cross- 
section (b) for 239 Pu ) averaged over the Maxwell spectrum, on the Maxwell and Fermi spectra 
joining energy and 77 = 1.8 (see (7)). 




a) b) 

Figure 3: Temperature dependencies for fission cross-section (a) and radioactive capture cross- 
section (b) for 235 U : averaged over the Maxwell spectrum, on the Maxwell and Fermi spectra 
joining energy and 77 = 1.8 (see (7)). 




1000 2000 3000 4000 5000 6000 




1000 2000 3000 4000 5000 6000 



a) 



b) 



Figure 4: Temperature dependencies for fission cross-section (a) and radioactive capture cross- 
section (b) for g 38 £7, averaged over the combined Maxwell and Fermi spectra depending on 
Maxwell and Fermi spectra joining energy and 77 = 1.8 (see (7)). 
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of neutron fitting into the resonance area. Therefore according to the known expressions 
for 235 U determining the neutron reactions cross-sections behaviour depending on their 
energy E n for non-resonance areas we observe dependency for the averaged cross-sections 

The data on the averaged fission and capture cross-sections of 238 U presented at fig. 4 
on one hand show that the averaged fission cross-section for 238 U, due to a high fission 
threshold ^1 MeV (see fig. 4a), is almost insensitive to neutron spectrum hardening 
caused by the fuel temperature increase, and on the other hand - confirm the dependency 
of capture cross-section on temperature, because its resonance area is located as low as for 
239 Pu. Obviously in this case fuel enrichment with 235 £/ essentially makes no difference 
because the averaged cross-sections for 235 U, as described above, behave in a standard 
way. 

And finally we performed a computer estimate of thermal source density dependence 
qrp(r, $,T, t) (6) on temperature for different compositions of uranium-plutonium fissile 
medium with a constant neutron flux density, presented at fig. 5. We used the dependen- 
cies presented above at fig. 2-4 for these calculations. Let us note that our preliminary 
calculations were made not taking into account the change in composition and density of 
the fissile uranium-plutonium medium that is a direct consequence of constant neutron 
flux assumption. 

The necessity of such assumption is caused by the following. It is obvious that for 
reasonable description of neutron source density q^ (r, <3>, T, t) (6) temperature dependency 
a system of three equations must be solved, two of which correspond to neutron kinetics 
equation (flux and fluence) and system of equations for kinetics of nuclear density of 
parental and child nuclides (e.g. see [19, 20]), while the third corresponds to a heat 
transfer equation of (1) type. However some serious difficulties arise here, associated with 
the computational capabilities available . And here is why. 

One of the principal physical peculiarities of TWR is the fact [21] that fluctuation 
residuals of plutonium (or 233 U in Th-U cycle) over its critical concentration burnout 
for time comparable with reactor lifetime of neutron r n (x,t) (not considering delayed 
neutrons) or, at least comparable with the reactor period 1 T(x, t) (considering delayed 
neutrons), while new plutonium (or 233 U in Th-U cycle) is formed in a few days (or a 
month) and not at once . This means [21], that numerical calculation must be made with 
temporal step around 10 -6 -10 -7 in case of not taking into account the delayed neutrons 
and ^10 _1 -10° otherwise. At first glance taking into account the delayed neutrons ac- 
cording to [21] really "saves the day", however it is not always true. If the heat transfer 
equation contains a significantly non-linear source, then in the case of a blow-up mode, 
the temperature at some conditions may grow extremely fast and in 10-20 steps (with 
time step 10 -6 -10 -7 s) reaches critical amplitude that may lead to (at least) a solution 
stability loss or (as maximum) blow-up bifurcation of the phase state, almost unnoticeable 
with rough time step. 

Here we should add that scientific literature has absolutely no theoretical or experimen- 
tal data on heat capacity c p and heat conductivity K of fissile material for temperatures 
over 1500K, which makes the further model calculations of the problem quiet problematic. 

According to these remarks, and considering the goal and format of this paper, we 
didn't aim at finding an exact solution of some specific system of three joined equations 
described above. However, we found it important to illustrate - at qualitative level - 
the consequences of possible blow-up modes in case of non-linear heat source existence 
in heat transfer equation. As said above, we made estimate computer calculations of 
thermal source density q^ (r, $,T, t) (6) temperature dependence in 300-1400K range for 
some compositions of uranium-plutonium fissile medium at constant neutron flux (fig. 5) . 

The obtained dependencies for source thermal flux density q J T (r, T, t) were success- 

x The reactor period by definition equals to T(x, t) = r n (x, t)/p(x, t), i.e. is a ratio of reactor neutron lifetime 
to reactivity. 
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Figure 5: Dependency of the thermal source density q* T (r, $, T, iV$, t) [eV] on the fissile medium 
temperature (300-1400K) for several compositions of uranium-plutonium medium (1 - 10% Pu; 
2 - 5% Pu; 3 - 1% Pu) at constant neutron flux density $ = 10 13 n/(cm 2 • s). 



fully approximated by a power function of temperature with an exponent of 4 (fig. 5). In 
other words, we obtained an equation for heat transfer with a significantly nonlinear heat 
source in the following form: 



where 5 > 1 in case of non-linear thermal conductivity dependency on temperature [4- 
8]. The latter means that the solutions of the heat transfer equation (1) describe the 
so-called Kurdyumov blow-up modes [4-9], i.e. such dynamic modes when one of the 
modeled values (e.g. temperature) turns into infinity for a finite time interval. As noted 
before, in reality instead of reaching infinite values a relevant phase transition is observed 
( final phase of parabolic temperature growth), that requires a separate model and is a 
basis for an entirely new problem. 

Mathematical modeling of blow-up modes was performed mainly using Mathematica 
5.2-6.0, Maple 10, Matlab 7.0, utilizing multiprocessor calculations for effective applica- 
tion. For calculations a Runge-Kutta method of 8-9 th order and the numerical methods 
of lines [22] were applied. The numerical error estimate didn't exceed 0.01%. The coordi- 
nate and temporal step were variable and chosen by the program in order to fit the given 
error at every calculation step. 

Below we give the solutions for heat transfer equation (1) with nonlinear exponen- 
tial heat source (11) in uranium-plutonium fissile medium for boundary and initial pa- 
rameters corresponding to technical reactors. The calculations were done for a cube 
of fissile material with different spatial size, boundary and initial temperature values. 
Since the temperature dependencies of thermal source densities were obtained without 
account for changing composition and density of uranium-plutonium fissile medium, dif- 
ferent blow-up modes can take place (HS-mode, S-mode, LS-mode) depending on a ratio 
between the exponents of temperature dependencies of thermal conductivity and heat 
source according to [4-9]. Therefore we considered cases for 1 st , 2 nd and 4 th temper- 
ature order sources. Here the power of the source also varied by varying the propor- 
tionality factor in (11) (const = 1.00 J/ (cm 3 • s • K) for 1 st temperature order source; 
0.10J/(cra 3 • s • K 2 ), 0.15 J '/(cm 3 • s • K 2 ) and 1.00 J /(cm 3 • s • K 2 ) for 2 nd temperature 
order source; 1.00 J /(cm 3 • s • K A ) for 4 th temperature order source). 

During calculations of the heat capacity c p (fig. 6a) and thermal conductivity K (fig. 6b) 



q T (T) = const -T^ 1+5) 



(11) 
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Figure 6: Temperature dependency of heat capacity cp and thermal conductivity x °f fissile 
material. Points - experimental values for heat capacity and thermal conductivity of 238 U. 



of fissile medium dependency on temperature in 300-1400K range the specified parame- 
ters were given by analytic expressions, obtained by approximation of experimental data 
f or 238 jj b asec [ on polynomial progression: 

c p (T) « -7.206 + 0.64T - 0.0047T 2 + 0.0000126T 3 + 

+ 2.004 • 10" 8 T 4 - 1.60 • 10" 10 T 5 - 2.15 • 10" 13 T 6 , (12) 

N(T) « 21.575 + 0.0152661T. (13) 

And finally the thermal transfer equation (1) solution was obtained for constant ther- 
mal conductivity (27.5 W/(m- K)) and heat capacity (11.5 J/(K -mol)) values, presented 
in fig. 7a, and also the solutions of heat transfer equation considering their temperature 
dependencies (fig. 7b-d). 

Preliminary results point directly to a possibility of local melting of uranium-plutonium 
fissile medium, having melting temperature almost identical to that of 238 U, that equals 
1400K (fig. 6a-d). Moreover, these regions of local melts are not the areas of so-called 
thermal peaks [23] and, probably, are the anomalous areas of uranium surface melt ob- 
served by Laptev and Ershler [24] that were also mentioned in [25]. More detailed analysis 
of probable temperature scenario associated with blow-up modes are discussed below. 

3 Physical peculiarities the blow-up regimes in 
neutron-multiplying media 

Earlier we noted the fact that due to coolant loss at nuclear reactors during Fukushima 
nuclear accident the fuel was melted or, in other words, temperature inside the active 
zone at some moment reached melting temperature of uranium-oxide fuel, i.e. ^3000°C. 

On the other hand, we already know that coolant loss may become a cause of non- 
linear heat source formation inside the nuclear fuel and therefore become a cause of 
temperature and neutron flux blow-up mode occurrence. A natural question arises, if it 
is possible to use such blow-up mode ( of temperature and neutron flux) for initiation of 
certain controlled physical conditions under which the nuclear burning wave would reg- 
ularly "experience" a so-called controlled blow-up regime. It is quite difficult to answer 
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Figure 7: Thermal transfer equation (1) solution for 3D case (crystal sizes 0.001x0.001x0.001 
mm; initial and boundary temperatures equal to 100K): a) The source is proportional to 4 th 
order of temperature; const =1.00 Jj (cm 3 • s • K A ) ) heat capacity and thermal conductivity 
are constant and equal to 11.5 J/(K • mol) and 27.5 W/(m • K) respectively; b) the source is 
proportional to 4 th order of temperature; const =1.00 J / (cm 3 • s • K 4 ); c) The source is propor- 
tional to 2 nd order of temperature; const =1.00 J/ (cm 3 • s • K 2 ); d) the source is proportional 
to 2 nd order of temperature; const = 0.10 J /(cm 3 • s • K 2 ). Note: in cases b)- d) heat capacity 
and thermal conductivity were determined by (12) and (13) respectively. 
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this question definitely, because such fast process has some physical vaguenesses, any of 
which can become experimentally unsur mount able for such process control. 

Nevertheless such process is very elegant and beautiful from the physics point of view 
and therefore requires more detailed phenomenological description. Let us try to make it 
in short. 

As we can see from plots of capture and fission cross-sections evolution for 239 Pu 
(fig. 2), the blow-up mode may develop actively at ^1000-2000K (depending on real value 
of Fermi and Maxwell spectra joining boundary), but at temperatures over 2500-3000K 
returns to almost initial cross-sections values. In case at that point we turn on some 
effective heat sink, the fuel may return to initial temperature state. However, while the 
blow-up mode develops, fast neutrons already penetrate to adjacent fuel areas, where new 
fissile material starts accumulating and so on. After some time in the this adjacent area 
the same blow-up mode will develop and everything starts over again. In other words, 
such hysteresis blow-up mode, closely time-conjugated to a heat takeoff procedure, will 
appear at background of stationary nuclear burning wave in a form of periodic impulse 
bursts. 

In order to demonstrate the marvelous power of such process, we investigated the heat 
transfer equation with non-linear exponential thermal source in uranium-plutonium fissile 
medium with boundary and initial parameters emulating the heat takeoff process. In 
other words, we investigated blow-up modes in fast Feoktistov-type uranium-plutonium 
reactor [2] where inside and outside boundary the temperature was deliberately fixed at 
6000K. 

This temperature is defined by an important question: "Is it possible in such conditions 
(6000K) emulating the time-conjugated heat takeoff (see fig. 2) to obtain a solution, i.e. 
spacio-temporal temperature distribution not in form of 5-function at some local spatial 
area, but as some kind of stationary and limited by amplitude solitary wave?" As shown 
below, such suggestion proved productive. 

Below we present some calculation characteristics and parameters. During these calcu- 
lations we used the following expression for dependency of thermal conductivity coefficient: 

K = 0.18-1CT 4 -T, 

which was obtained using Wiedemann-Franz law and data on electric conductivity of met- 
als at temperature 6000K [26] . Specific heat capacity at constant pressure was determined 
by value c p w 6 cal/(mol • deg) according to Dulong and Petit law. 

Fissile uranium-plutonium medium was modeled as a cube with size 10.0xl0.0xl0.0m 
(fig. 8) during calculations. Here for thermal source we used 2nd order temperature de- 
pendency (see (11)). 

And finally fig. 8a-d present a set of solutions of heat transfer equation (1) with non- 
linear exponential heat source (11) in uranium-plutonium fissile medium with boundary 
and initial conditions emulating such process of heat takeoff that initial and boundary 
temperatures remain constant and equal to 6000K. 

It is important to note here, that solution set, presented at fig. 8, is q demonstrate quite 
clearly the temporal evolution of the solution to its "Stationary" state. This is achieved 
using the so-called "magnifying glass" approach when the solutions of the same problem 
are deliberately investigated at different timescales. E.g. fig. 8a shows the solution at time 
scale t E [0, 10 -6 s], while fig. 8b describes spatial solution of the problem (temperature 
field) for t — 10 -6 s. The fig. 8c-d presents the solution - the spatial temperature 
distribution - at t = 0.5 s and t = 50 s. 

As one can see, the solution (fig. 8d) is completely identical to the previous (fig. 8c), 
i.e. to the distribution established in the medium in 0.5 seconds, that allowed us to make 
a conclusion on temperature field stability, starting from some moment. It is interesting 
that the established temperature field creates the conditions enough for thermonuclear 
synthesis reaction, i.e. reaching 10 8 K, and such temperature field existence time is not 
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Figure 8: Heat transfer equation solution for model reactor (source 2 nd order temperature 
dependency, const = 4.19 J/ (cm 3 • s • K 2 )] Initial and boundary temperatures equal to 6000K; 
fissile medium is a cube 10x10x10 m. The presented results correspond to the following times 
of temperature field evolution: (a) (1-10) -10 -7 s, (b) 10 -6 s, (c) 0.5 s, (d) 50 s. 
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less than 50 s. With a corresponding nuclei concentration to enter thermonuclear synthesis 
reaction, these condition, according to a known Lawson criterion, are favorable for a stable 
thermonuclear burning. 

However, one should remember that the results of this chapter are only demonstrative, 
because their accuracy is very relative and requires careful investigations with application 
of necessary computational resources. However, qualitative peculiarities of these solutions 
should attract the researchers' attention to the nontrivial properties of blow-up modes, at 
least, with respect to the obvious problem of internal TWR safety violation. 

4 Conclusion 

Below we give short conclusions stimulated by significant problems, that can be formulated 
in the following form. 

1. The consequences of anomalous 238 U and 239 Pu cross-sections behavior 
with temperature increase. It is shown that the capture and fission cross-sections 
of 238 U and 239 Pu manifest a monotonous growth in 1000-3000K range. Obviously 
such anomalous temperature dependency of 238 U and 239 Pu cross-sections cardinally 
change the neutron and thermal kinetics of nuclear reactors, TWR in particular. It 
becomes essential to know their influence on kinetics of heat transfer because it may 
become the cause of positive feedback with neutron kinetics, which may lead not 
only to undesirable loss of solution stability (the nuclear burning wave), but to a 
trivial reactor runaway with a subsequent nontrivial disaster. 

2. Blow-up modes and problem of nuclear burning wave stability. One of 

the causes of possible temperature growth of the fuel is, for instance, a deliberate 
or spontaneous coolant loss analogous to Fukushima nuclear accident. As shown 
above, the coolant loss may become a cause of nonlinear thermal source formation 
in nuclear fuel, and therefore emerging of mode with temperature and neutron flux 
blow-up . In our opinion, the preliminary results of heat transfer equation with non- 
linear thermal source investigations point to an extremely important phenomenon of 
anomalous behavior of thermal and neutron flux blow-up modes. This result poses 
a natural nontrivial problem of fundamental nuclear burning wave stability and, 
correspondingly, of a physically reasonable application of Lyapunov method to this 
problem, which is a basis for motion stability theory, and thus a reliable basis for 
justification of the Lyapunov functional minimum existence. 

It is shown that some variants of solution stability loss are caused by anomalous 
nuclear fuel temperature evolution. They can be not only the cause of TWR internal 
safety loss, but through a bifurcation of states (and this is very important!) can lead 
to a new stable mode when nuclear burning wave would periodically "experience" 
the so-called controlled blow-up regime. At the same time, it is noted that such fast 
(blow-up regime) process has a number of physical vaguenesses, any of which may 
happen to be experimentally insurmountable for control of such process. 

3. On-line remote neutrino diagnostics of intra-reactor processes. Due to the 

fact that a high-power TWR or a nuclear fuel transmutation reactor are projects 
of single-load and fuel burn- up with a subsequent burial of the reactor apparatus, 
there is an obvious necessity for remote system of neutrino control of neutron-nuclear 
burning wave in normal operation mode and control of neutron kinetics in emergency 
situation. The details and peculiarities of calculation of spatial-temporal distribution 
of isotope composition of active zone in TWR are presented in [27-29] in detail within 
the inverse problem of neutrino diagnostics of intra-reactor processes. 
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